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ABSTRACT 

We study the occurrence of delayed SNe la in the single degenerate (SD) scenario. We assume that 
a massive carbon-oxygen (CO) white dwarf (WD) accretes matter coming from a companion star, 
making it to spin at the critical rate. We assume uniform rotation due to magnetic field coupling. 
The carbon ignition mass for non-rotating WDs is MA R « 1.38 M 0 ; while for the case of uniformly 
rotating WDs it is a few percent larger (M R ps 1.43 M 0 ). When accretion rate decreases, the WD 
begins to lose angular momentum, shrinks, and spins up; however, it does not overflow its critical 
rotation rate, avoiding mass shedding. Thus, angular momentum losses can lead the CO WD interior 
to compression and carbon ignition, which would induce an SN la. The delay, largely due to the 
angular momentum losses timescale, may be large enough to allow the companion star to evolve to a 
He WD, becoming undetectable at the moment of explosion. This scenario supports the occurrence 
of delayed SNe la if the final CO WD mass is 1.38 M 0 < M < 1.43 M 0 . We also find that if the 
delay is longer than ~ 3 Gyr, the WD would become too cold to explode, rather undergoing collapse. 
Subject headings: binaries: close — nuclear reactions, nucleosynthesis, abundances — supernovae: 
general — stars: rotation — white dwarfs 


1. INTRODUCTION 

The origin of Type la supernovae (SNe la) has been 
one of the central issues of recent astrophysics. The main 
features of SNe la have been firmly accounted for as a 
consequence of a thermonuclear explosion of a carbon- 
oxygen (C+O) white dwarf (WD). Both Chandrasekhar 
mass and sub- Chan d rase kh ar mas s mod els have been 
prese nted ( Hillcbr andt fe Niemeverl 120001 : iNomoto et al.l 
120001) . However, it is still not clear if the the WD ac¬ 
cretes H/He-rich matter from its binary companion (sin¬ 
gle degenerate (SD) scenario), or t wo C+ O WDs me rge 
(double degenerate (DD) scenario) rtMaoz et al.ll2014 ). 

Observations provided some constraints on the nature 
of companion stars that may be considered as an indi¬ 
cation of the occurrence of the SD scenario. T hese are 
the presence of circumstellar matter 1CSM) (iPatat et al.1 

I2007t [Sternberg et al.ll201lllFolev et al .11 20121) and the de¬ 
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tection of hydrogen in the circu mstellar-interaction- type 
SNe (Ia/II n), such as SN 2002 ic (IHamuv et al .11 2 0031) and 
PTFllkx (IDildav et al.l 120121) . On the other hand, it 
has been not possible to detect the presence of com¬ 
panions, e.g., the lack of companion stars in the im¬ 
ages of S N 2011fe (iLi et al.ll201l[), so me SN la remnants 
(SNRs) dSchaefer &; Pamiottal 120 1 3) . SN 1572 (Tycho 
Kerzendorf et al.l 120091) . and SN 1006 (|Kerzendorf et al.l 
20121 ) . 

However, the above discussion does not take into ac¬ 
count that the accreting WD should be rotating by get¬ 
ting t he angul a r mo mentu m from the accre t ing m ateri- 
als. iJusthaml (|2011l) and iDi Stefano et al.l (120111) pre¬ 
sented a spin-down scenario and suggested that there 
should exist a long spin-down phase of the rapidly ro¬ 
tating WD with a timescale of the angular momen¬ 
tum loss (J-loss) from the WD. If so, the donor star 
in the SD model may exhaust its hydrogen-rich enve¬ 
lope and become a He W D before the SN la explo- 
sion. IHachisu et al.1 (|2012l) and IWang et~akl (|2014l ) stud¬ 
ied close binary evolutions by taking into account the 
rotation of accreting WDs and obtained the parameter 
space (the companion’s mass and the binary separation) 
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where the “ J-loss” evolution actually occurs. The struc¬ 
ture and evolution of rot ating WDs ha ve been studied b y 
many groups in ID (e.g.. lOstriker fc Bodenheimerll 19681 


L anger et al.l 2000) Piers anti et al 

2003; Saio & Nomotol 

20041: Yoon & Langeil 2001; IPiro 

120081) and 2D (Ag., 

Hachisull 19861 Uenishi et al.ll2003l 

Yoon & Laneedl2005 

)■ 


In the present study, we have constructed a new ID 
stellar evolution code by taking into account rotation. 
In Section [2] we describe the ingredients of the models 
as well as the numerical techniques employed. With this 
code, we have calculated the final evolution of the mass- 
accreting WDs by assuming that it rotates uniformly be¬ 
cause of magnet ic fiel d couplin g. By combining the close 
binary scenario (|Hachisu et al .1120111) , we have found that 
the WD evolves as follows, as it will be described in detail 
in Section [3] 

(1) For certain ranges of binary parameters, as it will 
be presented in Section [IJ the accretion rate (M) always 
exceeds 1O~ 7 M 0 yr -1 so that the WD increases its mass 
until it undergoes “prompt” carbon ignition. The mass of 
the uniformly rotating WD at the carbon ignition, M^, 

is larger for smaller M. For M = 10~ 7 M 0 yr -1 , Mj R = 
1.43 Mq, which is the largest mass because nova-like 
hydrogen flashes prevent the the WD mass from growing 
for the lower M. Because of the centrifugal force in the 
rotating WD, = 1 .43 M 0 is larger than M^ R = 
1.38 Mq dNomoto et alJlT984l) . 

(2) For adjacent ranges of binary parameters (JH), the 
mass of the rotating WD exceeds M^ g R = 1.38 M 0 but 
does not reach M : R = 1.43 M 0 because of the decreas¬ 
ing accretion rate. In this Letter, we shall present the 
evolution of the WD whose mass reaches 1.40 Mq, being 
slightly below M R g . 

After the accretion rate falls off, the WD undergoes the 
J-loss evolution. The exact mechanism and the timescale 
of the J-loss are highly uncertain, although we consider 
that the magneto-dipole braking WD to be responsible. 
We have found that J-loss induces the contraction of the 
WD, which leads to the “delayed” carbon ignition after 
the “delay” time due to neutrino and radiative cooling. 

We also studied how the carbon ignition depends on 
the strong screening factor, which still involves some un¬ 
certainties. Then, we shall discuss whether the WD un¬ 
dergoes the SN la explosion or collapse (jj3|). Finally, we 
discuss our findings in Section |4] 


2. METHODS 


We have empl oyed the hydros t atic s tellar evolution 
code described in lBenvenuto et ah (2013 ), adapted to the 
case of accreting WDs. The equation o f state o f degen¬ 
erate electrons is treated as in iKippenhahn &; Thomasl 
(11964f) . while Coulomb corrections at finite temperature 
are included following Hansc MMl Radiat ive opac¬ 
ities are taken from llglesias fc Rogersl (11996T) whereas 
condu ctive opacities were included following [Itoh et, al.1 
(2008) . Neut r ino em ission processes are those described 
in lltoh et al.l (il996ih 

As stated above, in this Letter we consider uniform 
rotation. For each isobar, we computed the parameter 
A, defined as 


Q 2 rp 

2 GM P ’ 


(1) 


where f 1 is the angular rotation, rp is the radius of a 
sphere with a volume equal to the corresponding isobar, 
and Mp is the embraced mass. A measures the fastness of 
rotation; breakup rotation rate corresponds to A « 0.270. 

Let us define the average < q > of an attribute q on a 
given isobar as 


<q>= 




( 2 ) 


where Sp is the isobar surface and da an element of it. 

In order to compute the structure of the rotating star, 
we need the coefficients fp and fp that correct the stan¬ 
dard structure equations for non-rotating objects. These 
are 


/ 47T7p A 1 

’ = \GMpSp) <gj>' 

47rr|> 

~s7 

where g e s is the effective gravitational accel erati on, com- 
puted in the Roche model approximation dMaederll2009H . 

When the WD is accreting, we impose the critical rota¬ 
tion regime for A = 0.270 for the outermost layers. This 
provides the (uniform) rotation rate of the entire star. 

For the mass accretion rate, we assume that M is 
constant during the main accretion stage and then de¬ 
clines exponentially with an assumed timescale r A ; f = 
4 x 10 4 years (set in order to mimic binary evolution) 
since a moment chosen in order for the WD to reach an 
assumed final mass value M. Since M < lx 10~ 7 M 0 
yr -1 , we also consider J-losses with another timescale 
tj, so that 



1 


< Jeff' > < Jeff > 


(3) 


J = J(t 0 )exp(-(t-t 0 )/Tj), 


(4) 


where to is the age at which M = lx 10 -7 M 0 yr -1 . 
We have chosen this prescription to explore the reac¬ 
tion of the WD to J-losses without assuming a particular 
physical process, which at present is uncertain. For this 
reason, we shall treat tj as a free parameter. A possible 
mec hanism leading to su c h losses is magneto -dipole brak¬ 
ing (lllkov fe Soker l2012t iRueda et al.lDOl.'f) . For t > t 0 
we compute the rotation rate by considering the instan¬ 
taneous moment of inertia of the oblate star. Due to 
J-losses, A falls down, although the WD may even rotate 
faster, leading to internal compression. 

For the carbon-burning 12 C+ 12 C reaction rate, we 
consi der the expressions given in iCaughlan &; Fowlerl 
(|1988f) . A key ingredient is the screening of this reaction. 
Here, we sh all explore four d if ferent t reatm e nts, which 
are th ose of lltoh et al.l (11990H. iQgatal (11997fl . 

(200C), and iPotekhin fe Chabried (120121) . 


Kitamura 


Kitamura 


1 20001 ) computed the enhancement factors of nuclear re¬ 
action rates due to interionic correlations within the 
linear-mixing (LM) approximation, where polarization 
of elec trons was explici tly taken into account. More re¬ 
cently, IPotekhin fc Chabrier ( 201 2) claimed that the de¬ 
viation Ah from the LM approximation tends to reduce 
the enhancement factor significantly, but their formula 
(19) is wrong; following their notations, the correct for¬ 
mula should read 
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Ah = h 0 - h\ m = 


/mix(^C; 2-0; 0) 


rf/mix(ac-(2-a;c)g,a;o(l+g),g) 

d€ 


«-*-0 


(5) 


where fmix(xc, ^Mg) represents the deviation of the 
interaction free energy from the LM value for C-O-Mg 
three-component plasma, with Xt (i= C, O, Mg) denot¬ 
ing their molar fractions. When this correct formula is 
adopted, Ah turns out negligibly small under the present 
parameter conditions. For practical use, they also sug¬ 
gested a simple Salpeter-Van Horn-type interpolation 
formula combined with the quantum correction, but we 
shall not use it here since it underestimates the reaction 
rates by one or two orders of magnitude near the freezing 
point, exhibiting large discrep ancies with the pycnonu- 
clear rat es in the so l id ph ase (IKitamural 120001) . We also 
note that lltoh et all (119901) and lOeatal (119971) ignored the 
electron polarization ef fects. Con s equen tly, we shall con¬ 
sider the treatment by IKitamural (I2000D as the standard 
case and is the one employed unless stated otherwise. 


3. RESULTS 

In this work, we restrict ourselves to the case of M = 
10 -6 M 0 yr -1 ; such material is burned in the outermost 
WD layers to C+O composition. The initial composition 
is ofX( 12 C) = X( 16 0) = 0.49 and heavier elements have 
solar abundances. 

We computed several evolutionary situations that de¬ 
pend on the binary parameters (orbital period P - donor 
mass M- 2 ) as shown in stopping calculations before the 
onset of hydrodynamical stages of thermonuclear run¬ 
away (log(T c /K) ~ 8.8: lNomotolll982j ). 

(1) One is the case in which the WD accretes material 
continuously up to “prompt” carbon ignition. As de¬ 
scribed in !JTJ we obtained = 1.43 Mq for M = 
10 -7 Mq yr -1 . This is smaller th an t h e carbon igni - 
tion mass of 1.48 Mq obtained by lUenislii et al.1 (j2003l ) 
because their 2D uniformly rotating model neglects the 
Coulomb term in the equation of state. 

(2) We also consider a binary system in which the ac¬ 
cretion rate declines and the WD reaches a final mass 
M before carbon ignition; in this work, we shall consider 
the case of M = 1.40 Mq. For t > to we apply J-losses 
as given in Equation (4) with timescales of tj = 1, 3, 
10, 30, 100, 300, and 1000 Myr. As quoted above, J-loss 
processes are highly uncertain; thus, we employ a wide 
range for tj to find the dependence of WD evolution 
with this parameter. The main results are presented in 
Figures fiTTl 

In Figure [1] we show the evolution of the spin period 
P sp in for the presupernova evolution of a WD. Accretion 
occurs for P sp m A) 2.9 s. After the end of accretion, the 
WD loses angular momentum exponentially (Equation 
(4)). Despite that, its moment of inertia decreases fast 
enough for the WD to spin up rather than spin down. 
We thus do not call such an evolution as “spin down,” 
rather just “J-loss” evolution. However, remarkably, the 
WD does not overflow the critical rotation rate avoiding 
shedding mass. This is shown in Figure [I] for the case 
of tj = 30 Myr. This is crucial for the delayed explo¬ 
sion scenario to work because a little mass loss would be 
enough to preclude the occurrence of the explosion. In 



Age [Myr] 


Fig. 1.— Spin period for the presupernova evolution o f a WD for 
the ca se of screening of carbon-burning reaction given by IKitamural 
(|2000I V The final mass is of 1.40 Mq. Accretion occurs for 
Pspin >2.9 sec. Each curve is labeled with its corresponding value 
of tj given in million years. Although the star loses angular mo¬ 
mentum exponentially, its moment of inertia decreases fast enough 
to make the star to spin up after the end of accretion. However, 
remarkably, the star does not overflow its critical rotation rate 
avoiding shedding mass. This is shown in the inset for the case 
of tj = 30 Myr, where spin and critical spin period evolution are 
given as solid and dashed lines, respectively. 


all cases, P sp i„ decreases (i.e., spin up rather than spin 
down) up to a minimum value slightly before explosion; 
however, A decreases monotonically. 

In Figure [21 we present the evolution of the center of 
the WD since before the end of accretion up to the on¬ 
set of the hydrodynamical stage. Thick black and thin 
red li nes corr espond to th e stro ng screeni ng fac tors g iven 
bv lKitamural (| 20001) and lPotekhin fe Chabrierl (120121) . re¬ 
spectively. The uppermost line corresponds to the case 
(1) evolution that leads to the “prompt” carbon igni¬ 
tion. Below that, the lines from upper to lower corre¬ 
spond to the evolutions for tj = 1, 3, 10, 30, 100, 300, 
and 1000 Myr, respectively, that lead to the “delayed” 
carbon ignition. 

In Figure [21 we also show the conditions at which neu¬ 
trino losses equals nuclear energy release (£ nuc = £u) and 
those that indicate when the central entropy per baryon 
begins to increase S = 0, i.e., carbon ignition. These 
curves are well detached because entropy per baryon 
starts to increases only when nuclear energy releases ex¬ 
ceeds neutrino and conductive losses. Conductive losses 
are comparatively more relevant the longer the tj be¬ 
cause neutrino losses are more sensitive to temperature 
than conductive losses. 
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Fig. 2.— Evolutionary tracks of the center of the WD up to the 
onset of the hydrodynamical explosion. The upper track corre¬ 
sponds to “prompt” carbon ignition, while the others are ordered 
from top to bottom with increasing J-loss timescale (tj = 1, 3, 
10, 30, 100, 300, and 1000 Myr, respectively). e n uc = £v indi¬ 
cates the conditions at which neutrino losses equal nuclear energy 
release, while S = 0 shows the stages at which central entropy 
per baryon begins to increase. Thick black and thin red lines cor ¬ 
res pond to the t reatments of sc reening given by [K itamura ( 2000) 
and Potekhin & Chabrier (2012), respectively. Arrows indicate the 
sense of the evolution. (A color version of this figure is available in 
the online journal.). 



Age [Myr] 


Fig. 3.— Evolution of the central temperature as a function of 
time for the cases included in Figure [T] 


In Figure [3] we depict the evolution of the central tem¬ 
perature as a function of time for the cases included in 
Figure [I] It is seen that for longer 7j, the WD is cooler 
and takes longer time until the carbon ignition occurs. 

In Figure [U we present the evolutionary tracks 
for the center of the WD from the carbon ignition 
for the cases of tj =3 and 100 Myr for the four 



Fig. 4.— Evolutionary tracks for the center of the WD during 
carbon ignition for the cases of tj = 3 and 100 Myr. I, O, K, 
and PC correspond t o the resul t s emp l oying the tre a tments for 
the sc reeni ng given in lltoh et al. (1990), Ogata (1997), Kit amural 
(2000), and [Pot ekhin & Chabrier (2012), respectively. Arrows in¬ 
dicate the sense of the evolution. 


strong scr e ening factors (I, O K, and PC) giv en in 
Itoh et al.j ( 19901 ), .Og ata) 019971) . iKitamural (l2000h . and 


Potekhin fe Chabried ( 20121) . respectively. Note that ig¬ 


nition densities are not so different among the four fac¬ 
tors. 

If tj > 3 Gyr, the central temperature, and thus the 
rate of nuclear energy release, becomes so low that car¬ 
bon is depleted before reaching the dynamical stage. 
Due to electron capture, such a case could lead to 
collapse of the WD rather than explosion, although a 
detailed study o f this phenomenon is necessary (e.g., 
iForster et al.lfeofoft . 


4. CONCLUSIONS AND DISCUSSION 

The calculations presented above show that a C+O 
WD accreting enough material coming from a binary 
companion, rotating uniformly at the critical rate, and 
reaching a mass M fulfilling M^ R < M < M R (i.e., 
1.38 Mq < M < 1.43M©) may indeed undergo delayed 
explosion. The above delay time is largely determined 
by the timescale of J-losses tj (see Figure [3]). For the 
values of tj considered here, the WD spends a time to 
undergo SN la explosion enough for the donor star to 
evolve to a structure completely different from the one 
it had when acted as a donor. For the red giant donor, 
its H-rich envelope would be lost as a result of H-shell 
burning and mass loss so that it would become a He WD 
in ~ 10 Myr. For the main-sequence donor, it would 






























































Single Degenerate Type la Supernovae 


5 



log P (day) 

Fig. 5.— Outcome of the binary evolution of the WD + compan¬ 
ion star systems is shown in the parameter space of the initial or¬ 
bital period P and the companion mass M 2 for the initial WD mass 
of 0.9Mq. The mass M of the WD starting from the “painted” re¬ 
gion reaches 1.38 Mq < M < 1.43 Mq (delayed carbon ignition), 
while the systems starting from the blank region encircled by the 
solid line reach M = 1.43 Mq (prompt carbon ignition). 


also evolve to become a low-mass He WD in 1 Myr, a 
hot He WD in 10 Myr, and a cold He WD in 1000 Myr 
(|Di Stefano et all 1201 ill . So, the J-losses should delay 
the explosion in enough time for the former donor to be 
undetectable. Therefore, this scenario provides a way to 
account for the failure in detecting companions to SNe la. 

Then, it is important to know for what binary systems 
the uniformly rotating WD effectively stops increasing 
its mass M at 1.38 Mg < M < 1.43 Mg and under¬ 
goes the delayed carbon ignition. We have calculated 
the binary evolution of the WD + companion star sys¬ 
tems with the parameters of the initial orbital period P 
and the companion m ass M 2 for the several initial WD 
masses as in lHachisu et all 2012f ). The results for the 
initial WD mass of 0.9 Mg is shown in Figure [5] Here, 
the binary systems starting from the “painted” region of 
the (P — M 2 ) plane reach 1.38 Mg < M < 1.43 Mg, 
while the systems starting from the blank region encir¬ 
cled by the solid line reach M = 1.43 Mg. We estimate 
that the occurrence frequency of the delayed carbon igni¬ 
tion would roughly be one-third of the total frequency of 
the carbon ignition. Th e occu r renc e frequ ency is smaller 
than the estimate by I Wang et all (l2014l) . possibly be¬ 
cause they adopted the mass range of the J-loss evolution 
as 1.38 Mg < M < 1.5 Mg. 

While here we have restricted ourselves to present re¬ 
sults corresponding to the case of WDs that reach a mass 
of 1.40 Mq with M = 10 -6 Mg yr -1 , we have explored 
this scenario for other final masses and accretion rates. 
We found that the overall evolution of the objects is qual¬ 
itatively the same if the final WD mass is close enough to 
Mh. While we did not attempt to determine the lowest- 
mass value for ignition to occur, we found that models 


with 1.37 Mg do not ignite, while those with 1.38 Mg 
explode. This is consistent with the adopted value of 

m| r . 

We explored different screening theories for the carbon¬ 
burning reaction finding that there is an uncertainty in 
the ignition density of « 10%. Remarkably, ignition (de¬ 
fined as the moment at which entropy per baryon of mat¬ 
ter at the center of the star begins to increase) occurs at 
a density significantly higher than the one at which nu¬ 
clear energy release equals neutrino losses. This is due 
to conductive losses, especially important for the case of 
large tj values for which the WD interior is cool. 

Carbon burning forms a convective core, which leads 
to the occurrence o f the convective URCA process (e.g., 
iLesaffre et all I2005T) . In the present study, we do not 
take into account this process, assuming that it does 
not suppress the heating by carbon burning. After the 
thermonuclear runaway, the central region would reach 
nuclear statistical equilibrium (NSE) and undergo elec¬ 
tron capture at high er den sities than in the non-rotating 
WD dNomoto et al.1 fTOSlll as well as the prompt igni¬ 
tion. This would lead to the synthesis of more neutron- 
rich Fe-peak elements, such as 54 Fe, 56 Fe, and 58 Ni. It 
would be important to study such nucleosynthesis in 
multi-dimensional simulation of convective deflagration, 
because the results would be sensitive to the treatment 
of turbulent flame. 

The nucleosynthesis results should be compared with 
observational signatures of SNe la. For example, the 
emission lines observed in the spectra of the nebular 
phase of SNe la indi cate t h e existen ce of a large amount 
of stable Fe and Ni. iMazzali et al.l (l2007f l estimated the 
amount of stable Fe is ~ 0.2Mg for most of normal SNe 
la (see also IMazzali et all 120151 for SN 2011fe). If such 
an amount of stable Fe is mostly 54 Fe and those SNe 
la are powered by the radioactive decay of ~ O.6M 0 
56 Ni, the ratio of 54 Fe/ 56 Fe in the ejecta of those SNe 
la is ~ 8 times larger than the solar ratio. This suggests 
that the stable Fe is mostly 56 Fe rather than 54 Fe, which 
would be realized for lower Y e , i.e., at the higher cen¬ 
tral dens ity than t he p rompt ignition model (such as W7 
in lNomoto et al.l I1984T) . Lower Y e in the delayed carbon 
burning might also be consistent with the large amount 
of Mn and (stable) Ni ar e observed in the SN la remnants 
(|Yamaguchi et al1l2015l) . 
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